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ABSTRACT 

Explicit  solutions  to  the  Cauchy  problem  for  the  linearized  KdV 
equation  are  constructed  when  the  initial  data  is  integrable.  The  method 
is  analogous  to  the  Fourier  decomposition  for  a  constant  coefficient 
equation  and  uses  the  connection  between  the  one-dimensional  Schrodinger 
equation  and  the  KdV  equation,  as  discovered  by  Gardner,  Greene,  Kruskal , 

and  Miura  [2] .  An  expansion  theorem  expressing  any  integrable  function  in 
terms  of  derivatives  of  squared  Schrodinger  (generalized)  eigenfunctions 

is  proved.  These  functions  evolve  according  to  the  linearized  KdV  equation, 
hence  the  expansion  of  the  initial  data  leads  to  a  generalized  solution 
of  the  linearized  KdV  equation.  Under  suitable  restrictions  on  the  initial 
data,  the  solution  constructed  is  classical.  The  proof  of  the  expansion 
theorem  may  be  interpreted  as  the  skew-adjoint  analogue  of  the  more  familiar 
process  of  simultaneously  diagonalizing  two  self-adjoint  operators. 
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SIGNIFICANCE  AND  EXPLANATION 


The  Korteweg-deVries  equation  (KdV  for  short)  arises  as  an  approxi¬ 
mation  in  many  non-linear  wave  problems  with  weak  dispersion  and  weak  non¬ 
linearity.  We  present  a  method  for  constructing  explicit  solutions  of  the 
linearized  KdV  equation.  This  equation  is  of  importance  in  studying  the 
effects  of  perturbations  of  the  KdV  equation.  The  method  relies  on  the 
connection  between  the  KdV  equation  and  the  one-dimensional  SchrcSdinger 
equation,  as  discovered  by  Gardner,  Greene,  Kruskal,  and  Miura  [2).  An 
expansion  of  the  initial  data  in  terms  of  derivatives  of  squared  (general¬ 
ized)  eigenfunctions  of  the  SchrBdinger  equation  provides  a  decomposition 
of  the  solution  resembling  the  use  of  Fourier  transforms  in  solving  constant 
coefficient  equations.  For  the  linearization  about  the  zero  solution  of 
the  KdV  equation,  the  analogy  holds  in  the  strict  sense,  as  our  expansion 
reduces  to  the  usual  Fourier  transform. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


COMPLETENESS  OF  DERIVATIVES  OF  SQUARED  SCHRODINGER 
EIGENFUNCTIONS  AND  EXPLICIT  SOLUT’IONS  OF  THE  LINEARIZED  KdV  EQUATION 


Robert  L.  Sachs 


1.  INTRODUCTION 

In  this  paper,  we  present  an  explicit  solution  to  the  Cauchy  problem 
for  the  linearized  KdV  equation: 


(*) 


u  +  u 

t  XXX 


-  6(qu)x 


0 


u(x,0)  =  <j>  (x) 


where  q(x,t)  is  a  solution  of  the  KdV  equation  (  (1.5)  below).  Our 
method  expresses  the  solution  as  a  superposition  of  particular  solutions 
and  utilizes  a  completeness  theorem  which  we  discuss  below.  The  particular 
solutions  we  choose  may  be  thought  of  as  derivatives  of  q(x,t)  with 
respect  to  the  scattering  data  for  the  Schrodinger  equation  with  potential 
q  (x , t) .  Hence  we  sketch  briefly  the  inverse  scattering  method  of  solving 
the  KdV  equation,  as  discovered  by  Gardner,  Greene,  Kruskal,  and  Miura  [2). 
If  we  consider  the  one-dimensional  Schrodinger  equation  with  potential 

Q(x)  , 

d2  o 

(1.1)  - J  f  +  Q(x)f  =  K  f 

dx 

and  define  the  Jost  solutions  f+(x,k)  by  their  asymptotic  behavior 

, ,  _ ,  ,  ikx 

(1.2)  f  e  as  x  *  +,x* 

,  -  ikx 

r  e  as  x  ►  -  ■■ 

then  the  relation 

(1.0  T(k)  f_(x,k)  =  ff(x,-k)  +  R (k )  f  (x,k) 
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which  defines  the  transmission  coefficient  T(k)  and  the  reflection  coef¬ 
ficient  R(k)  implies  T(k)  is  meromorphic  in  Im  k  >  0  with  finitely 

many  poles,  all  on  the  imaginary  axis.  The  completeness  theorem  mentioned 

2  2 

above  expresses  any  integrable  function  $  in  terms  of  (f  )  '  (x,k) ,  (f_  )'(x,k) 

and  a  sum  of  discrete  terms  related  to  the  poles  of  T(k).  (While  we  could 

2 

use  (1.3)  to  eliminate  (f_  ) • (x,k) ,  it  is  more  convenient  not  to  do  so.) 

We  prove  the  theorem  by  solving  the  equation 

tp1 '  ’  -  4Q4>'  -  2Q'  i>  +  4k 2\p'  -  4> 

for  and  integrating  the  'resolvent*. 

If  we  now  consider  a  one-parameter  family  of  Schrddinger  operators 
d2 

(1.4)  L(t)  = - r  +  q(x,t) 

dx 

where  the  time  evolution  of  q(x,t)  is  given  by  the  KdV  equation 
(1*2)  qt  +  qxxx  -  6qqx  =  0 

then  it  turns  out  [2]  that  for  any  t,  L(t)  is  unitarily  equivalent  to 
L (0) .  This  implies  that  the  spectrum  of  L(t)  is  invariant.  Moreover, 
the  scattering  data  associated  with  the  operators  L(t)  evolves  in  a 
very  simple  manner.  Zakharov  and  Faddeev  [15]  interpret  the  above  facts 
in  the  context  of  completely  integrable  Hamiltonian  systems  and  show  that 
the  eigenvalues  of  L  and  k  times  the  logarithm  of  |T(k)  |  for  k 
real  form  action  variables  with  appropriate  conjugate  angle  variables. 

In  [15] ,  a  formal  calculation  appears  which  expresses  the  infinitesmal 
variation  of  q(x,0)  in  terms  of  variations  of  the  scattering  data.  This 
formula  suggests  consideration  of  x-derivatives  of  the  squared  eigen¬ 
functions  of  (1.1)  with  their  induced  time  dependence  as  solutions  to  the 
linearized  KdV  equation.  The  fact  that  these  derivatives  satisfy  the 
linearized  KdV  equation  for  smooth  potentials  already  appears  implicitly 
in  [2],  Theorem  3,6.  A  related  expansion  for  the  Zakharov-Shabat  eigen- 
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value  problem  appears  in  Kaup  [4] .  Discussions  of  perturbations  using 
the  inverse  scattering  formalism  appear  in  [5] ,  [9]  (for  the  sine-Gordon 
equation) ,  and  [10] ;  for  an  application  of  this  result  to  the  problem  of 
water  waves  in  a  canal,  see  [12] .  Squared  eigenfunctions  and  their  deriv¬ 
atives  also  play  an  important  role  in  the  theory  of  the  periodic  KdV 
equation  [8] . 

The  completeness  theorem  is  proved  in  Section  2  below,  while  in 
Section  3,  the  time  evolution  of  the  eigenfunctions  and  the  solution  of 
(*)  are  discussed.  Some  of  these  results  appear  in  the  author's  doctoral 
dissertation  (N.Y.U.,  October  1980).  The  advice  and  encouragement  of  his 
advisor,  Jurgen  Moser,  is  gratefully  acknowledged. 

We  also  remark  that  semi-group  methods  [3]  will  yield  a  solution  of 
2 

(*)  for  L  initial  data,  so  tor  a  large  class  of  initial  data,  we  have 
constructed  the  "evolution  operator"  explicitly. 
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2.  L1-COMPETENESS  OF  DERIVATIVES  OF  SQUARED  SCHR&DINGER  EIGENFUNCTIONS 

After  introducing  some  notation  and  results  from  the  scattering  theory 
of  the  one-dimensional  SchrSdinger  equation,  we  state  and  prove  an  expan¬ 
sion  theorem  for  derivatives  of  squared  SchrSdinger  eigenfunctions.  (We 

shall  use  the  term  eigenfunction  to  include  generalized  eigenfunctions  as 

2 

well  as  bona  fide  L  solutions.) 

Consider  the  SchrSdinger  equation 

(2.1)  -  f”(x,k)  +  Q(x)  f(x,k)  =  k2f(x,k) 

for  k  real.  Our  notation  shall  be: 


f 1  (x,k)  =  f  (x,k) 


f  (x,k)  =  f (x,k) 


We  assume  the  potential  Q(x)  satisfies 


(2.2)  || Q  ||  =  /  (1  +  x  )  |Q(x)  |  dx  < 

L-  — 


The  fundamental  discovery  of  Gardner,  Greene,  Kruskal,  and  Miura  [2],  later 
formulated  abstractly  by  Lax  [7],  is  that  if  q(x,t)  evolves  according 
to  the  KdV  equation,  the  spectrum  of  the  Schrodinger  equation  (2.1)  with 
potential  q(x,t)  remains  fixed  in  t  and  the  associated  scattering  data 
evolves  in  a  very  simple  manner.  We  shall  use  this  information  below, 
but  first  introduce  some  notation  and  basic  facts  about  scattering  theory 
for  (2.1).  This  information  (and  much  more)  may  be  found  in  fl). 

Let  f+(x,k)  denote  the  Jost  solutions  of  (2.1) — 

—  ikx 

i.e.  f+(x,k)  ^  e  as  x  -+  +°°,  f  (x,k)  ^  e  as  x  * 

and  both  satisfy  (2.1).  The  transmission  coefficient,  T(k),  as  defined 

in  (1.3)  above,  is  represented  in  terms  of  the  Wronskian  of  f  ,  f  by: 


(2.3) 


f'f  -  f'f 

=  tf  (x,k)  ,  f  (x,k))  =  -  +  "  ' 


T(k)  2ik  1  + 


2ik 
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Formula  (2.3)  and  the  normalization  of  f  ,  f  imply  that  T(k)  is 

meromorphic  in  the  upper  half-plane  Im  k  >  0  with  poles  at  k  =  i£ ^ , 

2 

j  =  1,.  .  .  ,  N  where  each  energy  -8^  is  a  bound  state  energy  in  (2.1). 

N  is  finite  by  a  classical  estimate  assuming  (1+ | x | ) | Q (x) |  is  integrable 
T(k)  is  also  continuous  and  nonzero  for  real  k  ^  0.  For  notational 
ease,  we  also  introduce  for  j  =  1,  .  .  .  ,  N  the  following  pair  of  func¬ 
tions  : 

(2.4)  F.  (x)  =  f+2(x,ig.)  ;  G  (x)  =  c  f+ (x, i$ . ) -g . (x) 


.  ,  F  f  (X'iM  1 

-  T-  f  (x,k)  -  f  (x,k) 

l  dk  I  -  f  (x,i£.)  + 

+  3  - 


where  g^ (x) 


k=i8  . 
3 


and  c.  is  chosen  so  that  f  F.'(x)  G.(x)  dx  =  1  for  j=l,.  .  . ,  N.  The 
j  3  3 

expansion  theorem  mentioned  above  is: 

Theorem  2.1  Suppose  Q(x)  satisfies  (2.2).  If  <t>  (x)  is  continuous  and 
in  L1 ,  then 

°°+ie  ?  °° 

(2.5)  (a)  <fi  (x)  =  lim  /  T  (k)  .  /  K(x,y,k)  My)  dy 

e-K3-<*>+ie  -co 

N  « 

+  J  /  [F!  (x)  G  .  (y)  -  G!  (x)  F  .  (y)  ]  My)  dy 
j=l  -®  3  3  3  3 

and 


2 


(2.5)  (b)  <t>  (x)  =  /  /  t  (f^)’(x,k)  (y,k)  -  (ff ) »  (x.k)  (y , 


k)  I  :  fy ) 


—  oc  »co 


N  03 


+  l  /  [F  .  '  (x)  G  .  (y)  -  G  .  '  (x)  F  .  (y)  ]  My)  dy 

j  =  l  3  3  3  3 


where  the  kernel  K(x,y,k)  is  defined  as: 


h(x,k)h(y ,k) )  for  y  <  x 

f2 (x,k) f2 (y ,k)  for  y  >  x 

with  h(x,k)  f  f+(x,k)f  (x,k). 

We  remark  that  for  q  =  0,  f+(x,k)  =  e  and  the  above  expansion 

reduces  to  the  ordinary  Fourier  transform.  The  latter  representation  (2.5)  (b) 

is  more  convienient  in  applications  while  (2.5)  (a)  is  central  to  the  proof 

of  the  theorem.  Before  presenting  the  proof,  we  discuss  the  choice  of  the 

particular  kernel  K(x,y,k)  of  (2.6). 

If  f  and  g  are  both  C2  solutions  of  the  Schrodinger  equation  (2.1) 

2 

for  the  same  energy  k  ,  then  their  product  f  •  g  is  a  solution  of  the 
third-order  equation: 

(2.7)  <1)'"  -  4QtJd  -  2Q’\l>  =  -  4k  V 

Two  linearly  independent  solutions  of  (2.1)  generate  three  independent 

2  2 

solutions  of  (2.7)  e.g.  f  ,  fg,  g  .  We  choose  f  (x,k)  ,  T(k)  •  f  (x,k). 

Then,  solving  the  inhomogeneous  form  of  (2.7)  for  a  function  <J>  (x)  by 
variation  of  parameters  leads  to  an  expression  for  ij.>  in  terms  of  <j> . 

The  kernel  K(x,y,k)  is  the  'Green's  function'  for  this  problem.  Differ¬ 
entiating,  we  obtain  formally 

(2.8)  ip’  =  (D2  -  4(Q  -  k2)  -  20'D-1)"1^ 

which  we  integrate  as  though  it  were  a  bona  fide  resolvent  and  obtain  a 
multiple  of  the  identity.  The  calculation  of  this  integral  is  the  content 
of  the  following: 

Lemma  2.2.  Let  be  the  semicircle  in  the  upper  half-plane  of  radius 

R  traversed  from  -R  to  R.  Thi'n: 
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(2.9) 


lim 

R-XX' 


2ni  j 


T2(k) 


K(x,y,k)  <J)  (y)  dy 


}  dk 


=  4>(x)  for  all  4  which  are  continuous  with  4'  ^  L  . 

Postponing  the  proof  of  Lemma  2.2  momentarily/  we  show  that  Lemma  2.2 
implies  Theorem  2.1. 

Proof  of  Theorem  2.1,  given  Lemma  2.2:  Apply  Cauchy's  theorem.  For 

2 

Im  k  >  0  ,  the  integrand  in  k  has  poles  only  at  the  poles  of  T  (k) ; 

an  easy  calculation  shows  that  for  k  =  i8,  »  .  .  .  ,  iB.,  (recall,  the 

IN 

2  2  2 

bound  state  energies  are  -8.,  <  -8.,  -8,  <  0)  ,  the  pole  of 

N  N-l  1 

T(k)  is  simple  [1].  Thus  the  integrand  we  consider  has  double  poles  at 
k  =  ifj . ,  j  =  1 ,  .  .  . ,  N. 

For  k  =  i8.  ,  there  is  a  constant  a_.  such  that  f  (x,i8J  = 
tjf+(x,i6 J .  Thus  the  quantities: 


(2.10) 


f+2(x,k)  f_2(y,k)  -  h(x,k)  h(y,k)  and 


f+(x,k)  f_(y,k)  -  f_ (x,k) f+ (y,k) 


vanish  identically  in  x,y  when  k  =  i6j. 


Let  A.  be  the  residue  of  T(k)  at  k  =  i8..  Then  the  residue  of  the 
J  3 

left-hand  side  of  (2.9)  at  k  =  i8^  is  precisely: 


(2.11) 


2 n i  iB-  dk 
J 


K(x,y,k)  <}>(y)  dy 


k=i8. 


A2 

2iri  6  . 
J 


_0_ 

dx 


f+ (x,i8j)  •  f+ (y , i6 j ) g j 


(y) 


-  f+  (x,i6 .)g. (x)  f+(y,i6^) 


4(y)  dy 
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_ld 

where  we  recall  g.(x)  =  —  —  [f  (x,k) -a . f  (x,k) ] 

1  i  die  “  n  ■  k  1  p .  • 

3 

The  right-hand  side  of  (2.11)  comes  from  replacing  f  (x,iBJ  by 

crf+(x,i8_.)  after  using  the  remark  in  (2.10).  Now  we  deform  the  semicircle 

to  the  real  k-axis.  By  our  definitions  above,  the  deformation  contributes 
N  r°° 

the  terms  £  [F! (x)G  (y)  -  G!(x)F,(y)J  4>(y)  dy  which  appear  in  the 

j=l  3  3  31 

conclusion  of  Theorem  2.1,  in  addition  to  the  integration  along  the  real 
k-axis.  We  also  remark  that  despite  appearances,  neither  real  k  integral 
has  a  singularity  at  k=0.  This  follows  from  the  fact  that  either 
(i)  f+(x,0)  and  f  (x,0)  are  linearly  dependent,  which  by  the  same 
remark  as  above  in  (2.10)  implies  K(x,y,k)  vanishes  at  least  linearly 
in  k  as  k  -*■  0  in  Im  k  >  0,  or  (ii)  f  (x,0)  and  f  (x,0)  are 
linearly  independent,  so  by  (2.3),  T(k)  =  ak  +  o(k)  as  k  ■+  0.  In 
either  case,  there  is  no  pole  at  k  =  0.  (See  [1)  for  a  further  discussion 
of  phenomena  at  k  =  0  in  scattering  theory.)  Thus  we  have  proved  (2.5)  (a) 
assuming  Lemma  2.2.  To  obtain  (2.5) (b),  we  remark  that  the  difference  be¬ 
tween  the  two  k-integrals  integrates  to  0. 

The  difference  between  (2.5)  (a)  -  (2.5) (b)  is  precisely: 


OO 


oo 


1 

_GO 


dk 

4rrik 


K(x,y,k) 


<j>(y)  dy 


3  r  1 

where  K(x,y,k)  =  sgn(x-y)g^-  (k) f+ (x,k) f_ (y ,k)  -  T (k) f_ (x ,k) f+ (y ,k) 

Using  (1.3)  to  eliminate  f_(x,k),  f_(y,k),  we  see  that  K(x,y,k)  is  an 
even  function  of  k  ,  so  the  integral  vanishes.  This  proves  (2.5) (b) . 


To  complete  the  argument,  we  now  prove  Lemma  2.2.  Consider 


(2.12) 


R  2iTi 


T2(k) 


K(x,y,k)  <My)  dy 


}  dk 
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We  shall  make  use  of  the  following  estimnt 


Write  f+(x,k)  =  m+(x,k)  e 

taken  from  Deif t-Trubowitz  [1]  ,  which  hold  for  all  k  in  Im  k  >  0 : 

C 

(2.13)  (i)  |m+(x,k)-l|  <  exp  {C^/ 1  k  |  }  y^y 

(ii)  |m|(x,k)  |  <  C3./(l+|k|  ) 

(iii)  T(k)  =  1  +  0(y£y)  as  M  “• 

In  [1]  it  is  also  shown  that  m+-l  are  Hardy  functions;  in  particular, 
they  are  analytic  in  Im  k  >  0.  Recalling  the  definition  of  K(x,y,k) 
given  by  (2.6),  we  define: 


(1) 


2ni  J 


T2(k)  / 

dk  — -  { 


i  1 


2  2 

2ik  m  (x#k)m  (y,k) 


.  e2iMx-y)  ;(v) 


+  2ikm2  (x,k)m2 (y,k)e  2lk (x  a(. 


<fj(y )  dy 


(2) 


1 

)ti 


dk 


T2(k) 


[2m  (k,x)m' (x,k)m2(y,k)e2lk(X  v) 
+  +  - 


h’  (x,k)  h  (y,k)  ]  <J>(y )  dy 


f  rn  < 


+  [h1  (x,k)h  (y,k)-2m  (x,  k)m_^  (x,k)m  (y,k)  •  e 


1  i(\ 


mu  r  T(l)  T (2) 

Thus  IR  =  IR  +  IR 


(2) 


By  estimates  (i) ,  (ii),  (iii)  in  (2.13),  since  I'  ’  contains  terms 

R 

which  have  a  factor  mj(x,k)  and  m  (x,k)  is  uniformly  bounded  for 
| k |  >  c  >  0  ,  we  have  the  estimate 


(2.15) 


U<2)!  <c- 

C  is  independent  of  R. 


for  all  R  _>  Rg  sufficiently  large,  where 


Moreover,  by  (2.13)  (i)  ,  (iii)  ,  m  (x,k)  =  1  +  O(~0.  T(k)  =  1  +  O(^) 

+  R  K 


for  k  =  R  so  I 


(i) 


i 

2-ni 


dk 

k 


2ik(x-y)  ,  ,  ,  , 

2ik  e  1  <P  (y)  dy 


...  -2ik(x-y)  ,  . 

2ik  e  4>(y)  dy 


*  0(r’ 

The  first  terms  converge  to  4>(>0  as  in  the  usual  proof  of  Fourier  complete¬ 
ness  [14]  and  the  lemma  is  proved  by  taking  R-*°°. 

Remarks :  (i)  The  expansion  theorem  above  bears  a  strong  resemblance  to  that 

of  the  Fourier  transform  for  L1  functions.  However,  since  the  underlying 

process  is  the  "simultaneous  diagonalization"  of  the  two  skew  operators 

— -  and  -  (-r~)  3  +  2—  Q+2Q-r~,  the  analoque  of  the  Fourier  l/  theory 
dx  dx  dx  dx 

does  not  exist  if  Q  ^  0.  If  we  define 


(2.16)  <f+(k)  = 


<My)f+(y,k)  dy  , 


the  natural  version  of  the  Plancheral  formula  in  this  case  relates  the 

oo 

skew  bilinear  form  f  i|t'(x)  <f>(x)  to  the  standard  symplectic  pairing 


/ 

^  (k> 

0  1 

<J>+  (k) 

(k) 

f 

1 

M 

O 

(k) 

(ii)  We  also  note  that  the  skew  operator  -D  +  2D Q  +  2QD  (D  =  — ) 

was  used  by  Lenard  to  recursively  generate  the  KdV  conservation  laws  [2] . 
It  seems  to  play  a  crucial  role  in  many  aspects  of  the  KdV  theory  —  e.g. 
it  is  useful  in  proving  that  the  integrals  are  in  involution. 
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3.  APPLICATION  OF  THEOREM  2.1  TO  THE  CAUCHY  PROBLEM  FOR  THE  LINEARIZED 


The  expansion  of  4>  (x)  given  by  Theorem  2.1  will  lead  directly  to 
a  method  for  solving  the  Cauchy  problem  for  the  linearized  KdV  equation: 

(*)  U  +  u  -  6  (qu)  =  0;  u(x,0)  =  <Mx) 

L  XXX  X 

with  <j)  satisfying  the  hypotheses  of  Theorem  2.1.  As  the  potential  q(x,t) 
in  the  SchrSdinger  equation 

-  f"  +  qf  =  k2f 

evolves  according  to  the  KdV  equation,  the  corresponding  eigenfunctions 
evolve  in  time.  Gardner,  Greene,  Kruskal,  and  Miura  [2]  observed  that 
the  squares  of  the  eigenfunctions  satisfy  the  formal  adjoint  of  the  linear¬ 
ized  KdV  equation  (which  they  called  "the  associated  linear  equation") , 
namely: 

(3.1)  v  +  v  -  6qv  =  0 

t  XXX  X 

from  which  it  follows  that  u  =  v  satisfies 

x 

(3.2)  u^  +  u  ~  6(qu)^  =  0  ,  the  linearized  KdV  equation. 

In  view  of  this  fact,  the  expansion  of  Theorem  2.1  may  be  extended  to 
include  the  time  evolution  of  the  eigenfunctions.  As  we  shall  see  below, 
this  extension  is  the  solution  of  the  Cauchy  problem  (*) . 

We  begin  by  developing  some  necessary  preliminary  facts: 

Lemma  3.1  The  functions  g.(x)  =  -7-  37-  [f  (x  ,k)  -a  .  f ,  (x,k)  ] ,  k  =  i£S.  are 
-  3  1  dk  3  +  3 

2  2 

(unbounded)  solutions  of  the  Schrbdinger  equation  (2.1)  with  k  =  -6^  . 
Proof :  Differentiating  (2.1)  with  respect  to  k  ,  we  obtain  an  equation 
for  (d/dk) f  (x ,k)  =  f+(x,k): 


(3.3) 


Consider 


(q  +  ej2)gj 

(q  +  Bj2)gj 


f;  =  (q  -  k2)f+  -  2kf+ 

-  a^f+]  at  k  =  igj.*  (x) 

-  28^  (f_(x,ig..)  -  a_.f+(x,i$  J  ) 
by  our  choice  of  a , . 


satisfies 


Remark:  g_.  (x)  is  exponentially  increasing  as  |x|  however,  the 

product  f+(x,i3j)  g^ (x)  is  bounded. 

We  now  discuss  the  result  in  [2]  mentioned  above,  and  sketch  the  proof. 
Lemma  3 . 2  (cf.  equation  (2.19)  of  [2]). 

Let  ip  be  a  solution  of  the  Schrbdinger  equation  (2.1)  with  potential  q(x,t) 
evolving  according  to  the  KdV  equation. 

Then  the  function 


(3.4) 


5  *t  + 


ip  -  3(q  +  k 

XXX 


is  also  a  solution  of  the  Schrbdinger  equation  with  potential  q(x,t). 

Sketch  of  the  proof:  Use  the  equation  ip"  =  (q-k2)^  to  express  q  in 

terms  of  i p.  Substitute  into  the  KdV  equation  and  simplify  the  resulting 

2 

expression,  obtaining  the  equation  R"  -  (q-k  )R  =  0  after  eliminating 
a  factor  of  ip . 

Remark :  The  chief  use  of  Lemma  3.2  is  to  show  that  R,  for  suitable 
eigenfunctions  ,  is  in  fact  0.  In  particular,  we  have: 

Corollary  3.3  The  expression  R  vanishes  if  we  choose  any  of  the  follow¬ 
ing  eigenfunctions  for  'i': 


(i)  f+(x,k,t)  'v  exp  j+i  (kx+4k^t) j  as  x  -*•  +  <»,  t  fixed; 

(ii)  f+(x,i0j,t)  ; 

(iii)  qAx,t)  :  t-  ~  (f_(x,k,t)  -  n  .  f+  (x,k ,  t ) ) k_.^ 
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Sketch  of  proof  (see  [2],  Theorem  3.6):  Consider  R  in  each  case.  1  / 
our  choice  of  asymptotics,  we  have  R  =  0  since  no  other  solution  of  the 
Schrbdinger  equation  with  that  type  of  decay  exists,  namely  R  -*■  0  as 
x  ->  4^  »  in  (i)  ,  exp{  |  8 ^x  -  48^t|}R-*-  0  in  (ii),  (iii)  as  1 8 jX-48  2t  j  >•  ■  . 
Thus  from  the  spectral  theory  of  the  Schrbdinger  equation,  R  =  0  in  each 
case. 

Lemma  3.4  Suppose  iji  ,iji  are  two  (not  necessarily  independent)  solutions 

of  the  Schrbdinger  equation  for  the  same  eigenvalue.  If  ^  +  'i>xxx 
2 

-  3  (q  +  k  =  0  for  tp=  j  =  1,2,  then  the  product  ^2^2  1S  3 
solution  of  the  adjoint  equation: 


v  +  v  -  6qv  =  0 
t  XXX  X 


Proof :  A  direct  calculation: 

(Wt +  (w*«  -  ^iVx 

=  ^ ^  ($2  +  "  6<3^2  )  +  ^2(’*'l  +^1  "  6qlyi  1 

t  xxx  x  t  xxx  x 

+  3<P1  i P2  +  3^ 

X  XX  XX  X 

=  4*-,^  (»^2  +  >1^  -  3(q+k2)^2  ) 

t  XXX  X 

+  ^  ($2  +  “  3(q+k2)(}<1  ) 

t  xxx  x 

=  0,  where  we  used  the  SchrSdinger  equation  to  eliminate 
the  second  derivatives. 

Remarks :  (i)  An  alternate  derivation  of  these  facts  may  be  given  using 

the  following  idea  of  Tanaka  [13] :  The  KdV  equation  may  be  written  in  the 
Lax  (7)  form 


(3.5) 


f-  IB'11 


where  L(t)  is  the  operator 


+  q(x,t)  and  B(t)  is  the  skew  operator 


dx 


(3.6) 


d  d  d 

_4  — _  +  3q(x,t'  gj  +  3  ^  q(x,t) 

dx 

The  time  derivative  of  the  Schrbdinger  equation  (2.1)  and  (3.5)  implies 

L(ffc  -  Bf)  =  k2(ffc  -  Bf). 

Choosing  f  =  f  (x,k,t)  'v  e  as  x  -*■  +  °°  for  t  fixed  and  analyzing  the 
asymptotic  behavior  of  f^  -  Bf  as  x  ->■  +  00  for  t  fixed  implies 

(3.7)  (f+)t  -  B(f+)  =  4 (ik) 3 (f+) 

so  ip(x,k,t)  =  e  •  f+(x,k,t) 

satisfies 

iVy  -+•  4j  V 

(3.8)  -  BiJj  =  0  ,  tj)  ^  e  as  x  -*-  +  °°  for  t  fixed 

3 

-4ik  t 

(a  similar  argument  holds  for  e  *  f_(x,k,t)).  Using  the  Schrbdinger 

equation  (2.1),  a  simple  calculation  as  in  Lemma  3.4  above  shows  that 

2 

products  (with  the  same  values  of  k  )  of  solutions  of  (3.8)  satisfy  the 
adjoint  equation 


v^  +  v  -  6qv 

t  XXX  X 


0. 


(ii)  If  q(x,t)  is  a  classical  solution  of  the  KdV  equation,  the  formal 
calculations  above  are  sensible —  i.e.  the  eigenfunctions  possess  the 
necessary  derivatives.  This  follows  from  the  inhomogeneous  form  of  the 
Schrbdinger  equation  which  these  derivatives  satisfy. 

Using  these  facts,  we  make  the  following  definitions  (extending  (2.4) 
to  include  time  dependence) : 
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(3.9) 


t+(x,k,t)  ^  e 


ikx+4ik3t 


as  x  -*■  +  60 ,  t 


~  ,  ,  . .  -ikx-4ik  t 

f  (x,k,t)  e  as  x  t 


f  satisfy  (2.1)  with  potential  q(x,t) 


(x,t)  =  f+(x,i6j,t) 


G_.  (x, t)  =  Cj  f+(x,iB  ,t)  •  g..(x,t) 


fixed 


fixed 


The  obvious  candidate  for  the  solution  of  (*)  is  the  function  u(x,t) 
defined  as  follows: 


(3.10) 


00 

u(x,t)  =  |  d-k^.(kk)  [(fj)  *(x,k,t)  $_(k)  -  <fV(x,k,t)  $+(k)] 

— <X> 


+ 


N 

I 


j=i 


CO 

[F  '  (x,t)G  (y,0) 

— oo 


Gj  ’  (x,t)F_.  (y , 0)  ]  <{>  (y)  dy 


where  (k )  are  defined  in  (2.16)  above.  By  Lemma  3.4,  all  of  the 
functions  of  (x,t)  appearing  on  the  right-hand  side  of  (3.10)  satisfy 
the  linearized  KdV  equation.  Thus  we  have  proved: 

Lemma  3 . 5  The  function  u(x,t)  defined  by  (3.10)  above  satisfies  the 
linearized  KdV  equation: 


(*) 


+  u 

XXX 


-  6(qu)x 


0 


u  (x,o)  =  <j>  (x) 


in  the  sense  of  distributions. 

To  prove  that  u(x,t)  is  a  classical  solution  of  (*)  for  t  ?  0  , 
we  need  some  additional  smoothness  and  decay  on  (x)  ,  the  initial  data. 
The  situation  is  completely  analogous  with  the  Fourier  transform  solution 
of  the  Airy  equation: 
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(3.11) 


V7  +  W  =0 
t  xxx 


The  x-decay  of  4>  is  needed  to  have  $+(k)  be  differentiable  while  smooth¬ 
ness  of  4>  relates  to  integrability  of  ka$^(k)  for  0  <  a  <  2.  As  in  the 
case  of  (3.11),  one  can  prove  the  following: 

Theorem  3 . 6  The  function  u(x,t) ,  given  by  (3.10)  above,  is  a  classical 
solution  of  the  linearized  KdV  equation  for  t  >  0  if 

(i)  <j>(x)  has  four  continuous  derivatives. 

(ii)  As  |x|  -*  »,  3^<j>(x)  =  0 ( ]  x |  4)  for  r  =  0,1,2, 3,4. 

Sketch  of  proof;  Using  (i)  and  the  definition  of  $+(k),  we  integrate 
by  parts  four  times  with  respect  to  x  ,  which  implies,  as  |k|  -*  00 , 

I  •  —A  A  2 

<fi+(k)  =  0(|k|  ).  Also,  by  (ii) ,  4>  +  (k)  are  C  .  Thus  u(x,t)  has  two 

continuous  derivatives  with  respect  to  x.  As  in  Murray  [11] ,  the  factor 

2 

2  +  12k  t  x 

k  is  written  as  - - ^2t  '  where  the  first  term  is  a  multiple 

— 2ik (x+4k^t) 

of  the  k-derivative  of  the  exponentials  e  .  Integrating  by 

parts  in  k  ,  we  find  that  u(x,t)  has  four  continuous  x-derivatives. 
Repeating  the  argument,  we  obtain  six  x-derivatives.  From  the  linearized 
KdV  equation,  this  implies  ufc  is  continuous,  hence  it  is  a  classical 
solution.  The  only  difference  between  this  case  and  the  Airy  equation 
(3.11)  is  the  presence  of  the  added  factors 

~  ,  ...  ~  .  ...  +2ik(x+4k2t) 

m  (x ,k , t)  =  f+(x,k,t)  e 

and  these  do  net  affect  the  necessary  estimates. 

A  fuller  discussion  and  an  analysis  of  asymptotic  behavior  is  presented 
for  the  N-soliton  linearization  in  [12],  where  the  perturbation  theory  for 
the  problem  of  water  waves  in  a  canal  is  discussed.  The  KdV  equation  was 
first  derived  to  model  precisely  this  situation  [6]. 
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